Self-sustained magnetoelectric oscillations in magnetic resonant tunneling structures 
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The dynamic interplay of transport, electrostatic, and magnetic effects in the resonant tunneling 
through ferromagnetic quantum wells is theoretically investigated. It is shown that the carrier- 
mediated magnetic order in the ferromagnetic region not only induces, but also takes part in intrinsic, 
robust, and sustainable high-frequency current oscillations over a large window of nominally steady 
bias voltages. This phenomenon could spawn a new class of quantum electronic devices based on 
ferromagnetic semiconductors. 



Ferromagnetism of diluted magnetic semiconductors 
(DMSs) , such as GaMnAs Jl| , depends strongly on the 
carrier density 0, 0, 0, H, The possibihty to tailor 
space charges in semiconductors by bias or gate fields 
naturally suggests similar tailoring of magnetic proper- 
ties of DMSs. While early experiments have indeed suc- 
ceeded in generating ferromagnetism in DMSs electrically 
or optically 0, S], ramifications of the strong carrier- 
mediated ferromagnetism in the transport through DMS 
heterostructures are largely unexplored. 

In resonant tunneling through a quantum well not only 
the tunneling current, but also the carrier density in the 
well are sensitive to the alignment of the electronic spec- 
tra in the leads and in the well. If the quantum well 
is a paramagnetic DMS, the resulting transport is influ- 
enced by the spin splitting of the carrier bands in the 
well, as observed experimentally 0]. The magnetic res- 
onant diodes are prominent spintronic device s l id . 
proposed for spin valves and spin filtering 13, 13l|. or for 
digital magnetoresistance 
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If the quantum well is made of a ferromagnetic DMS 
0, 121, resonant tunneling conditions should influence 
magnetic ordering as well. It has already been pre- 
dicted that the critical temperature Tr of the well can 
be strongly modified electrically fisl. Il9l. [20| . Here we 
show that the magnetic ordering affects back the tun- 
neling current, in a peculiar feedback process, leading to 
interesting dynamic transport phenomena. 

Conventional nonmagnetic RTDs can exhibit subtle 
intrinsic bistability and terahertz current oscillations 



21|, 122, 1231, |2J, |25| resulting from the nonlinear feedback 
of the stored charge in the quantum well. Interesting phe- 
nomena occur also in multiple quantum wells and super- 
lattices, in which electric field domains form whose dy- 
namics leads to current oscillations in the kHz-GHz range 
[i^t . This effect has been exploited for spin-dependent 
transport by incorporating paramagnetic quantum wells 



27, 2 



In this article we introduce a realistic model of a 
self-consistently coupled transport, charge, onrf magnetic 
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FIG. 1: (Color online) (a) Schematic scheme of the band pro- 
file of the magnetic double barrier structure. The exchange 
interaction of the magnetic ions is mediated by the carriers 
tunneling in and out of the well. Here, the cut off of the emit- 
ter tunneling rate is realized by a cascaded left barrier, (b) 
Equivalent circuit model of the resonant tunneling structure 
introducing the emitter and collector capacitances Ce, Cc and 
resistances R^,Rc, respectively. 



dynamics and apply it to generic asymmetric resonant 
diodes with a ferromagnetic quantum well to predict self- 
sustained, stable high-frequency oscillations of the elec- 
tric current and quantum well magnetization. We for- 
mulate a qualitative explanation for the appearance of 
these magnetoelectric oscillations. In essence, ferromag- 
netic quantum wells exhibit strong nonlinear feedback 
to the electric transport since the ferromagnetic order, 
which gives rise to the exchange splitting of the quantum 
well subbands, is itself mediated by the itinerant carriers. 
This together with the Coulomb interaction which effec- 
tively modifies the single electron electrostatic potential 
in the well, leads to a strong coupling of the transport, 
electric, as well as magnetic properties of ferromagnetic 
resonant tunneling structures. 

Our model resonant tunneling structure with a ferro- 
magnetic well made of a DMS-material, e.g., GaMnAs, is 
sketched in Fig.[TJa). To exhibit magnetoelectric oscilla- 
tions the structure needs a built-in energy cut-off, E'cut, 
of the emitter tunneling rate. Such an energy cut-off 
might be realized by a cascaded left barrier, as shown in 
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Fig. [TJa), by which the tunnehng for carriers with ener- 
gies smaller than Ecut is exponentially suppressed, due to 
the increased barrier width. Other possibilities to realize 
the cutoff are discussed below. 

As we aim to understand the most robust features of 
the ferromagnetic resonant tunnel structures, we present 
a minimal theoretical model which captures the essen- 
tial physics. The longitudinal transport through the sys- 
tem can be described in terms of sequential tunneling, 
as the high densities of magnetic impurities residing in 
the ferromagnetic well will likely cause decoherence of 
the propagating carries. Based on the transfer Hamil- 
tonian formalism a master equation for the semiclassical 
particle distribution in the well can be derived, as de- 
scribed elsewhere l3, 111- We assume that there is only 
a single resonant well level Eq in the energy range of 
interest, allowing us to write the rate equations for the 
spin-resolved, time dependent quantum well particle den- 
sities nc{t), {a =t, i= ±1/2) as. 



drier 



(1) 



Here, {e.c} denotes the energy-dependent tunneling rate 
from the emitter (e) and the collector (c), F = Fg + Fc 
denotes the total tunneling rate, Ts stands for the spin 
relaxation time in the well, no.a denotes the quasiequilib- 
rium particle spin density, and no-^{e,c} ^-re the densities 
of particles in the emitter and collector reservoir, respec- 
tively, having the resonant longitudinal energy E„. The 
spin-split resonant energies are 



E„ = Eo + U^,- (tA, 



(2) 



with Uw being the electrostatic well potential and A de- 
noting the subband exchange splitting. The total en- 
ergy of the particles i^tot is then given by the sum of 
the longitudinal energy Ea- and the in-plane kinetic en- 
ergy: Etot — Ea + Tin- The physical meaning of the right 
side of Eq. ^ is as follows: the first two terms are the 
gain terms, describing tunneling from the emitter and the 
collector into the well; the third term describes all loss 
processes due to the tunneling out of the well, and the 
last term models the spin relaxation in the well. Con- 
sidering the Fermi-Dirac distributions in the emitter and 
the collector, the particle densities 

= DakBTlnll + exp [{fi, - i?cr)/fcsr]} i = e, c, 

(3) 

with fcs denoting Boltzmanns' constant; T is the lead 
temperature, [ii are the emitter and collector chemical 
potentials with /ic — /^e— eT4ppi where T4ppi is the applied 
bias, and Dq = m/2TTh^ is the two-dimensional density 
of states per spin for carriers with the effective mass m. 
The tunneling rates are essentially given by the overlap of 



the lead and well wave functions according to Bardeen's 
formula [30| . For high barriers the rate becomes propor- 
tional to the longitudinal momentum pz of the particles 
[igt - i.e., Fe,c cx (i?^)^/^ with E^ denoting the longitudi- 
nal energy. 

In the framework of a mean field model for the carrier 
mediated ferromagnetism in heterostructure systems 
H, IE 113] the steady state exchange splitting of the well 
subbands is determined by 



(4) 



Ao = Jpd / dz nimp(z) |V'o(2)r 



xSBs 



SJpdsjni - n|) |^o(z)|' 



Here, Jpd denotes the coupling strength between the im- 
purity spin and the carrier spin density (in case of GaM- 
nAs p-like holes couple to the d-like impurity electrons), 
z is the longitudinal (growth) direction of the structure, 
n-impiz) is the impurity density profile, ipoiz) labels the 
well wave function, Bs denotes the Brillouin function 
of order S, and S and s = 1/2 are the impurity and 
particle spin, respectively. (In the case of Mn impuri- 
ties S — 5/2.) The expression shows that the well spin- 
splitting depends basically on the particle spin polariza- 
tion ^ = n| — n| in the well and the overlap between the 
wave function and the impurity band profile. For simplic- 
ity, we consider here a homogenous impurity distribution 
in the well, which makes the determining factor for Ao. 

After a sudden change of the well spin polarization 
the magnetic impurities need some time to respond until 
the corresponding mean field value Aq is established. In 
the case of GaMnAs, experimental studies of the mag- 
netization dynamics revealed typical response times of 
about 100 ps [3l[. We model the magnetization evolu- 
tion within the relaxation time approximation, dA/di — 
— (A — Ao)/rA, with ta denoting the well spin-splitting 
relaxation time. 

Finally, to take into account the nonlinear feedback of 
the Coulomb interaction of the well charges, we introduce 
emitter- well and collector- well capacitances, C = Ce + Cc 
according to the equivalent circuit model of a resonant 
tunneling diode, as shown in Fig.[l](b). The capacitances 
Ce and C'c are determined by the geometrical dimensions 
of the barriers and the well [2^. The electrostatic well 
potential can then be written as 



[e'^{n - nback) - C'ccVappi] 



(5) 



with e denoting the elementary charge, rtback is the pos- 
itive background charge (from magnetic donors) in the 
well. All the equations are nonlinearly coupled via Eq. ^ 
for the resonant well levels E^ , making a numerical solu- 
tion indispensable. 

In the numerical simulations we use generic parame- 
ters assuming a GaMnAs well: m = 0.5 niQ, Sr = 12.9, 
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FIG. 2: (Color online) (a) Current-voltage (IV)-characteristic 
of the investigated structure at the energy cut-off _Ecut = 85 
meV. In the unsteady region between 11 to 22 mV the current 
does not reach a steady state value; instead intrinsic current 
oscillations appear, in which both the maximum and mini- 
mum values of the oscillations are indicated by the dashed 
and solid lines, respectively. Two different dynamic modes I 
and II can be identified, (b) Well splitting A versus applied 
voltage. In the unsteady region the maximum (dashed) and 
minimum (solid) values of the magnetization oscillations are 
plotted, (c) and (d) Transients of the spin-resolved current j 
and the well particle density n at the applied voltage V = 12 
mV and -Ecut ~ 85 meV. The time is measured in units of t* 
(typically some picoseconds) as explained in the text. 



50 A, 4 = 20 A, w = 10 A, 
= 1.5 X 1020cm-3^ Jpd = 0.06 eV nm^, t, = It* 
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100 meV, 



TA = lot*, where de,dc and w are the emitter barrier, 
collector barrier and quantum well widths, toq denotes 
the free electron mass, and is the relative permittiv- 
ity of the well. The characteristic time scale t* is the 
inverse emitter tunneling rate at the emitters' Fermi en- 
ergy, [t* = l/VeifJ-e)], being of the order of picoseconds. 
For the well background charge we consider that in GaM- 
nAs the actual carrier density is only about 10% of the 
nominal Mn doping density [32|: JT-back = 0.1 riimp. Our 
calculations are performed at T — 4.2 K, which is well 
below the critical temperature of our specific well Tc ~ 10 
K, where we estimated Tc by using the mean-field result 
given in Ref. [20| . 

Figure [SJa) shows the current- voltage (IV) character- 
istic of the structure. Up to about V = 11 mV the 
typical peaked IV-curve of a resonant tunneling diode 
is obtained. However, in the voltage range of 11 to 22 
meV, which we will call hereafter the "unsteady" re- 
gion, the current does not settle down to a steady state 
value; instead stable high-frequency oscillations occur, 
as shown in Fig. ^c) for the applied voltage of ^ = 12 
mV. The current is always evaluated at the collector side: 
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lations of the well magnetization as well as of the spin 



FIG. 3: (Color online) Explanation for the occurrence of self- 
sustained oscillations in the case of mode I. For mode II the 
arguments are completely analogous but there the spin down 
level is crossing the cut-off energy -Ecut. (a) According to 
the in-tunneling spin-up carriers the spin polarization ^ is 
increased, giving rise to an increasing well splitting A. (b) 
When the spin up level falls below the cutoff energy the total 
particle number decreases and, hence, also the electrostatic 
potential Um does, (c) This pushes the spin-up level even 
deeper into the cut-off region leading to a fast decrease of the 
spin polarization and consequently of the well splitting, which 
brings the spin up level back to the emitter's supply region 
(d), restarting the whole cycle. 



densities appear, as shown in Fig [2jb) and (d). Those 
magneto-electric oscillations are the main results of this 
paper. 

The IV-curve in the unsteady region suggests the ex- 
istence of two qualitatively different dynamic modes. In- 
deed, comparing the transients in these two voltage re- 
gions reveals that in region I the spin up level is recur- 
ringly crossing the emitter energy cut off i^cut , whereas in 
region II this is done by the spin down level, as schemat- 
ically illustrated in Fig. [3l This insight offers the follow- 
ing explanation for the the occurrence of self-sustained 
oscillations. Take mode I; the arguments for mode II are 
similar. The dropping of the spin up level below the cut- 
off energy (due to the increasing exchange splitting) as 
two implications: (i) the supply of the emitter spin up 
electrons sharply decreases. Hence, the total well par- 
ticle density n = n| + rt; decreases because the spin up 
electrons residing in the well are tunneling out to the col- 
lector. A decreased particle density leads to a decreased 
electrostatic potential according to Eq. (O, which effec- 
tively drives the spin up level even deeper into the cut-off 
region, (ii) Since the spin up electrons are the majority 
spins in the well, a decrease of implies a decreasing 
spin polarization ^ in the well. This causes, via Eq. (j4|), a 
rapid decrease of the subband exchange splitting, bring- 
ing the spin up level back to the emitter supply region. 
The spin up electrons can then tunnel again into the well 
and the whole process starts from the beginning, produc- 
ing the calculated cycles. 

The occurrence of these oscillations needs the concur- 
rent interplay of both the electric and magnetic feedbacks: 
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FIG. 4: (Color online) Contour plots of the frequency (a) and 
the period (b) of the intrinsic oscillations as a function of 
the applied voltage and the cut-off energy -Bcut- The two is- 
lands of high-frequency oscillations correspond to the dynamic 
modes I and II. They are separated by a crossover region III 
of high-period oscillations, which becomes most evident in 
period contour plot (b). 



the electrostatic feedback acts like an "inertia" for the os- 
cillations, allowing the spin level to get deeper into the 
cut-off region, whereas the magnetic feedback is needed 
to bring the level back into the emitter's supply region. 
From the above discussions it also follows that a steep 
descent in the tunneling rate at the cut-off energy is nec- 
essary. This is confirmed in our simulations, where we 
assumed an exponential decay of Te for Ez < Ecut- Ac- 
tually, this is also the reason why no oscillations are oc- 
curring at the emitter's conduction band edge {E^ — 0), 
which is an otherwise natural emitter's cut-off, because 
near the band edge the tunneling rate is already almost 
vanishing according to Te oc (iJ^)^/^. 

The typical timescale for the oscillations is made up by 
two contributions: (i) the evolution of the well splitting 
A in the emitter's supply region, i.e., the time needed 
for the well splitting A to become large enough that one 
of the spin levels crosses the cut off energy, and (ii) the 
dynamics of A in the cut-off region. The timescale of 
contribution (i) can be estimated to be proportional to 
T^(r — 1), where r is the ratio of the well splitting when 
the cut-off energy is reached to the initial splitting at the 
beginning of the cycle, whereas the dynamics of contri- 
bution (ii) is mostly governed by a fast rearrangement 
of the spin densities in the well, which happens on the 
order of the tunneling time t*. This is in line with the 
numerical result that the oscillation frequency increases 
with decreasing ta. In contrast, a decreasing spin relax- 
ation time, which diminishes more and more effectively 
the spin polarization in the well and, hence, the well split- 
ting, gives rise to a decrease of the frequency. 

The fastest magnetoelectric oscillations in our simu- 
lations for the generic parameters used exhibit periods 
of about 2.5 x t* « 5 ps. Important, the oscillation fre- 
quency can be modified by varying the applied voltage, as 
can be seen in Fig. IH which displays contour plots of the 
frequency and the period versus the applied voltage and 
the cut-off energy. Two "islands" can be distinguished, 
corresponding to the two regions of the dynamic modes 



I and II, respectively (see Fig. [2]). Mode II appears at 
higher voltages, where the collector chemical potential 
is already below the cut-off energy. In this case the de- 
crease of the well particle density is much stronger when 
the spin up level crosses the cut-off energy, as illustrated 
in Fig. OJb) , since none of the out-tunneling electrons are 
"Pauli-blocked" by the collector electrons. This gives rise 
to a much stronger decrease of the electrostatic potential 
as compared to mode I, pushing both the spin up and 
down levels below the cut-off energy. In the following 
oscillations this makes the spin down level recurringly 
crossing the cut-off energy instead of the spin up level as 
is the case for mode I. 

Both modes are separated by a crossover region III, in 
which the oscillations have much longer periods than on 
the "islands", as can be seen in Fig.[4ljb). The reason for 
these low-frequency oscillations is that the initial spin po- 
larization and, hence, the well splitting at the beginning 
of each cycle is small, yielding a large ratio r. Therefore, 
in region III it can take 10-100 times longer that the well 
splitting reaches the cut-off energy as for modes I or II. 
The possibility of controlling the frequency of the mag- 
netoelectric oscillations by electrical means is especially 
interesting from the applications point of view. 

For the experimental observation of the self-sustained 
oscillations the emitter needs to emit electrons in a 
sharply defined energy window. In the model above we 
propose an asymmetric barrier, as in Fig. [Ija). Another 
practical way would be using an auxiliary resonant tun- 
neling structure in the emitter itself, providing a spectral 
filter allowing only resonant electrons to pass through. 
Yet another possibility would be employing a hot elec- 
tron emitter, say using a Schottky barrier. 

We have shown that the charge and magnetization 
dynamics in ferromagnetic tunneling heterostructurcs 
are influenced by the highly nonlinear feedback of both 
Coulomb and magnetic couplings on the tunneling trans- 
port. The feedback results in high-frequency self- 
sustained intrinsic oscillations suggesting applications of 
ferromagnetic QWs in tunable high-power current oscil- 
lators. 

This work has been supported by the Deutsche 
Forschungsgemeinschaft SFB 689. 
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